Theory of Raman response in a superconductor with extended s-wave symmetry: 

Application to Fe-pnictides 
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We argue that Raman study of Fe-pnictides is a way to unambiguously distinguish between various 
superconducting gaps proposed for these materials. We show that A\ g Raman intensity develops 
a true resonance peak below 2A if the pairing gap has A\ g symmetry in the folded Brillouin zone 
(A(fc = 0) = A, A(n, 7r) = —A). No such peak develops for a pure s-wave gap, a d-wave gap, and 
an extended s-wave gap with A(k) = A cos 4p- cos We show that the peak remains quite strong 
for the values of inter-pocket impurity scattering used to fit NMR data. 



PACS numbers: 74.20.Mn, 74.20.Rp, 74.25.Jb, 74.25. Gz 

Recent discovery of superconductivity in the iron- 
based layered pnictides with T c reaching 55K generated 
an enormous interest in the physics of these materials^. 
Most of fcrropnictides are quasi two-dimensional mate- 
rials, and their parent (undoped) compounds are met- 
als and display antifcrromagnctic long-range order be- 
low Tn ~ 150Ki2i£. Superconductivity occurs upon 
doping of cither electrons or holes into FeAs layers, or 
by applying pressure. The electronic structure mea- 
sured by angle-resolved photoemission (ARPES)^ and by 
magneto-oscillations^ consists of two small hole pockets 
centered around the T = (0, 0) point and two small elec- 
tron pockets centered around the M = (n, w) point in 
the folded Brillouin zone (BZ). The sizes of electron and 
hole pockets are about equal in parent compounds. 

The key unresolved issue for the pnictides is the 
symmetry of the superconducting gap. A conventional 
phonon-mcdiated s-wave superconductivity is unlikely 
because electron-phonon coupling calculated from first 
principles is quite small^. Several authors considered 
magnetically mediated pairing based either on itiner- 
ant 7 -^ or localized spin models^ and argued that the 
gap should have an extended s-wave symmetry cos k x + 
cos k y (also called s + or, equivalcntly, Ai g symmetry). 
This gap changes sign between hole and electron pock- 
ets but has no nodes along the Fermi surface (FS). On 
the other hand, another RPA study of magnetically medi- 
ated superconductivity in the five-band Hubbard model^i 
yielded two nearly degenerate candidate states in which 
the gap has nodes on one of the FS sheets: either an 



A cos -f cos or a 



extended s-wave state with A(k) 

d x 2_ y 2 state with A(k) tts A sin sin 4^ (in the unfolded 
BZ, these two states are cos q x + cos q y and cos q x — cos q y , 
respectively^) . 

The experimental situation is also controversial. 
ARPES 15 d 6 and Andreev spectroscopy^ 7 , measurements 
have been interpreted as evidence for a nodeless gap, ei- 
ther pure s-wave or s + -wave. The resonance observed in 



s+ garj±2. 



neutron measurements 1 ^ is consistent with the 
On the other hand, nuclear magnetic resonance (NMR) 
data— and some of penetration depth data3i were inter- 
preted as evidence for the nodes in the gap. Some, but 
not all of the data can be reasonably fitted by a model 
of an s + SC with ordinary impuritic o 14 i 22 i 23 . 

In view of both theoretical and experimental uncertain- 
ties, it is important to find measurements which could un- 
ambiguously distinguish between different pairing sym- 
metries. Recent suggestions for such probes include An- 
dreev bound stated and Josephson interferometry^. In 
this communication, we argue that the study of Ai g Ra- 
man intensity is another way to determine the symmetry 
of the superconducting gap. We show that in the A\ g 
scattering geometry the Raman signal develops a true 
resonance below 2A for the case of s + gap. No such res- 
onance appears for a pure s-wave gap, for cos 4p cos 

and sin 4f sin ^ gaps. The A lg resonance is the effect of 
the final state interaction, which is known to be impor- 
tant for Raman scattering^. A similar resonance occurs 
in the B\ g channel in a magnetically mediated d x i_ y n 
superconductor—, but there the resonance is weakened 
by a finite damping associated with nodes of the d-wave 
gap. 

We model Fe-pnictides by an itinerant electron system 
with two (almost) degenerate hole FS pockets centered at 
the r point and two electron FS pockets centered at the 
M point. We assume that the magnitude of the gap A 
is much smaller than the Fermi energy. In this situation, 
Raman intensity at frequencies < 2A is determined by 
states near the FS where the density of states (DOS) 
can be approximated by a constant. We first assume 
that the pairing gap has s + symmetry, A(k ps 0) = A, 
A(k ps 7r) = —A, and show how the resonance appears. 
We then discuss other pairing symmetries. 

Without final state interaction, the Raman intensity 
in a clean BCS s + superconductor is the same as in a 
pure s-wave superconductor— and is given by = 
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FIG. 1: Real and imaginary parts of the Ai g Raman intensity 
for a clean s + superconductor without (a) and with (b) final 
state interaction. Final state interaction gives rise to a well- 
defined resonance in the Ai g intensity. We used Ro = l/(47r), 
u £ ffRo « 0.4, and added damping 7 = 0.001 A. 



2lmRi(n), where 
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Here j Al = cos k x + cos fc y is the geometrical factor for 
A\ g scattering, uj± = uj ± 51/2, and (...) FS denotes the 
averaging over FS. The factor 2 in the relation between 
and Ri(£l) reflects the fact that there are two hole 
and two electron pockets. Other factors are incorporated 
into R . 

The intensity I Al (O) computed using ([I]) vanishes at 
£1 < 2 A and is discontinuous at 2 A. The real part of 
R Alg , which we will need later, is positive below 2A, 
scales as Q 2 at small frequencies, and diverges upon ap- 
proaching 2A from below^. We show both KeR Alg and 
Imi?^ in Fig. [T] The final state interaction is diagram- 
matically represented as the renormalization of the Ra- 
man vertex. Vertex corrections arise from multiple in- 
sertions of fermion-fermion interactions into the Raman 
bubble. There are five different interactions between low- 
energy fcrmions (see Fig.[2ja)). They include intra-band 
interactions for electrons (1/4) and for holes (us), which 
we assume to be equal, inter-band interactions u\ and «2 
with momentum transfer and (it, 7t), respectively, and 
the pair hopping term u 3 . 

A generic theory of vertex renormalizations has been 
developed in Ref3i and we follow this work in our anal- 
ysis. In general, there are three different types of ver- 
tex corrections: (i) the corrections which come from 
short-range interactions Ui and transform a bare A\ g 
Raman vertex into a renormalized particle-hole vertex 
(these corrections involve GG and FF bubbles, where G 
and F are normal and anomalous Green's functions), (ii) 
the corrections which transform a Raman vertex into a 
particle-particle vertex (these involve GF bubbles) (see 
Fig. [H D )), and (iii) the corrections from the long-range 
component of the Coulomb interaction V q cx 1/q 2 . The 
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FIG. 2: (a) Five relevant interactions between fermions near 
hole and electron FS pockets. Black and grey lines represent 
fermionic c- states near (0,0) and /- states near M = (ir,n). 
(b) The full Raman bubble, which is the sum of GG, FF and 
GF terms. Only the contribution from c— fermions is shown. 
l)The one from /-fermions is obtained by replacing c lines by 
/-lines and vice versa in (b) and (c) panel, (c) The renormal- 
ization of the Ai g Raman vertex for c- fermions. The first 
8 diagrams account for "conventional" renormalization of the 
Ai g particle-hole vertex and involve GG and FF bubbles, the 
last two diagrams involve GF bubbles and emerge due to a 
non-zero coupling between Ai g particle-hole channel and or- 
dinary s— wave pairing channel. The renormalization of the 
particle-particle vertex in turn involves 4 "conventional" dia- 
grams with GG and FF terms, which account for the renor- 
malization in the particle-particle channel, and two diagrams 
due to the coupling to the Ai g particle-hole channel, (d) The 
renormalizations due to long-range component of Coulomb 
interaction V q oc 1/q 2 . This renormalization vanishes because 
of the symmetry between c— and /—fermions and the fact 
that 7A lg (fc = 0) = -7-4 l3 W- 



corrections of the first type are given by ladder and bub- 
ble diagrams which involve u\, 112 and M4 vertices (first 
8 terms in Fig.f^c)). The corrections of the second-type 
are non-zero when the symmetry of the gap is the same 
as the symmetry of the Raman vertex, which is our case. 
These corrections transform A\ g particle-hole vertex into 
an ordinary s— wave pairing vertex (the terms with the 
overall factor 2 in Fig.[2](c)). The third renormalization, 
due to long-range component of the Coulomb interac- 
tion, generally gives rise to a screening of the Raman 
signa l 26 ' 29 , but vanishes in our case because of particle- 
hole symmetry and the fact that A\ g Raman vertex 7a 1sj 
changes sign between hole and electron pockets (see Fig. 
[2{d)). Note in this regard that V q oc 1/q 2 is not a part 
of RG transformation and depends only on a momentum 
transfer q, in distinction to the other two interactions 
with small momentum transfer, U4 and U\. The bare 
values of U4 and u\ may be identical, but the two flow 
in different directions under RG. Also note that we did 
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not include a momentum-independent term into "f Alg ■ If 
jAi g had such component, it would be screened by the 
long-range Coulomb interaction. 

Combining the renormalizations (i) and (ii) and eval- 
uating the diagrams, we obtain the full Raman intensity 
l{ u lg (Cl) = 2InxR^(0) with R f A "(Cl) in the form 

<>) = 

R Alg (n) (i + u f R PP m + 4u f R 2 mix (n) 

(1 - u eff R Alg (Cl))(l + u f R pp (Cl)) - Au g u f Rl lx (n) ' 

(2) 

where u e ff = 2u\ — it 2 — U4 is the effective vertex for the 
Raman renormalization in the A\ g particle-hole channel, 
u f = W3+U4, u g = U4—U2, Rpp(Q) oc log-Bir/f2 > is the 
polarization bubble in the s— wave particle-particle chan- 
nel, and R mix (CI) oc / d 2 kduj^ Alg G k ^ + nF k ^ couples A lg 
particle-hole channel and s-wave particle-particle chan- 
nel. At low frequencies, R m i X (Cl) oc CI. In this respect, 
the situation is similar to the case of a spin resonance 
in a d— wave superconductor, where S — 1 particle-hole 
channel couples to S = particle-particle channel 3 -^. 

Because s-wave channel is repulsive in our case (it3 + 
U4 > 0), there is no pole in R A u ^ l (Cl) coming from the 
particle-particle channel. Furthermore, R pp (Ct) logarith- 
mically diverges at Cl << Ef, and canceling this diver- 
gent term between the numerator and the denominator 
in we obtain 

/£"(fi)«2 (3) 

(1 - u eff RcR Alg ) + (u e fflmR Alg ) 

We see therefore that the coupling between Ai g particle- 
hole and s— wave particle-particle channels is irrelevant, 
and the full I A ™ (Cl) can be approximated by the expres- 
sion which only includes vertex corrections which pre- 
serve particle-hole structure of the Raman vertex. 

Our next observation is that for two-band structure, 
u e ff contains the terms u\ and u%, which do not con- 
tribute to the renormalization of the s + pairing vertex 
(the latter involves 113 and U4 terms 1 ^), i.e., in distinc- 
tion to one-band case2&, the renormalization of the A\ g 
Raman vertex and the renormalization of the s + pairing 
vertex (which has the same Ai g symmetry) are given by 
different combinations of the interactions u. 

Finally, we note that below 2A, lmR Al = while 
HcR Alg is positive and evolves between zero and in- 
finity when Cl changes between zero and 2A. Then, 
for positive u e //, the A\ g Raman intensity develops a 
^-functional resonance peak below 2 A, at a frequency 
where u e f fReR Alg = 1. For a d x 2_ y 2 superconductor 
the same effect leads to an excitonic resonance in a stag- 
gered spin susceptibility 3 ^, and to a pseudo-resonance in 
a Bi g Raman response^. 

The flow of the interactions between the bandwidth W 
and the Fermi energy Ep has been analyzed in the earlier 
RG studyi^, and the result is that u\ becomes the largest 



interaction at energies comparable to the Fermi energy, 
even if the intra-band Hubbard repulsion 114 is the largest 
term in the Hamiltonian. Specifically, in the RG flow u\ 
and U3 increase, u-iju\ flows to zero, and the U4 term 
decreases, such that u e f / = 2u\— 114— U2 becomes positive 
at energies below Ep, relevant to Raman scattering, and 
the Ai g Raman response develops a resonance below 2A. 
We emphasize that the physics which makes u e f / positive 
is the same physics that gives rise to an attraction in an 
extended s + -wave pairing channel. Indeed, the pairing 
interaction in s + channel, which is the combination U3 — 
M4, becomes positive under RG. 

For other proposed gap symmetries, the resonance does 
not develop, even if one neglects the screening by long- 
range Coulomb interaction. For an s-wave gap, there is 
no sign change between electron and hole FS, and the 
analog of U e ff in Eq. ([3]) is — 2iti — U4 + U2- This combi- 
nation is negative, so the resonance does not occur. For 
a gap that changes sign along either hole or electron FS, 
the largest contribution to I Al (CI) comes from the FS 
along which the gap is nodeless. Vertex renormalization 
for such term contains — U4 + (2ui — U2)x, where x ~ k F , 
and kp is a small radius of the FS along which the gap 
has nodes. In this situation, u\ term does not overcome 
U4, and the resonance does not occur. For d xy gap— with 
A oc sin/c^ sinfcy (in the folded BZ), all Ui terms in the 
vertex renormalization are reduced. Resonance may still 
occur, but the effective interaction now is small, 0(k F ), 
and the resonance is washed out by a small damping. 
This shows that the Ai g Raman resonance is a finger- 
print of an s + pairing. 

Finally, we consider how the resonance in s + supercon- 
ductor is affected by ordinary impurities. As in earlier 
work o 14 ' 23 , we introduce impurity potential Ui(q) with 
intra- and inter-pocket terms Ui(0) and Ui(jr), respec- 
tively and restrict with the Born approximation. This 
approximation (which requires Ui << Ep) may not work 
for Ui(0) (Ref.— ) but should be valid for Ui(ir) which is 
pair-breaking and is therefore very likely not larger than 
A << Ep. For our case, U(0) controls the functional 
form of ReR Alg , which still evolves between zero and 
infinity when CI changes between and 2A, while the 
broadening of the resonance is entirely due to Ui(ir). In 
this situation, Born approximation should be sufficient. 

The calculations are straightforward and we refrain 
from presenting the details. Intra-pocket impurity scat- 
tering does not affect the gap by Anderson's theorem, but 
Ui(n), which scatter fermions with +A and —A, is pair- 
breaking and affects the gap in the same way as magnetic 
impurities in an ordinary s-wave superconductor. We use 
b = 2Ui(ir)/A, where A is the order parameter as a mea- 
sure of the strength of pair-breaking impurity scattering. 

The results of the calculations are shown in Fig. [31 
where we plot Raman intensity in the presence of im- 
purities both without and with final state interaction. 
Comparing this figure with Fig. [T] we see that the reso- 
nance gets damped at a finite 6, and Raman intensity 
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FIG. 3: (color online) Calculated Raman intensity for an s + 
superconductor without (a) and with (b) final state interac- 
tion for various strength of the inter-band impurity scatter- 
ing. We use the same it e // as in Fig. [T]and for definiteness 
set Ui(0) = A. 

no longer shows two peaks. Still, the resonance con- 
tinue to determine the shape of Ia 1£I (Q)'- without final 
state interaction the peak broadens and shifts to larger 
frequencies Q > 2A upon increasing b, while when the 
final state interaction is included, the peak remains be- 
low 2A and shifts to a smaller frequency with increas- 
ing b. Notice that the resonance is still quite strong at 
b ~ 0.5 — 0.7, which was used to fit NMR and penetration 
depth dat a 14 ' 23 . In other words, it should be observable 



in Raman experiments if indeed the gap has an s + sym- 
metry. 

To conclude, in this paper we argued that Raman 
study of Fe-pnictides is a way to unambiguously distin- 
guish between various superconducting gaps proposed for 
these materials. We have shown that for an A\ g (s + ) gap 
A(k « 0) = A, A(k ps n) « -A, the A lg Raman inten- 
sity has a true resonance peak below 2 A. No such peak 
emerges for a pure s-wave gap, a d x 2_ y 2 gap, and an ex- 
tended s-wave gap with A(k) = A cos 4f cos The res- 
onance peak gets broader by pair-breaking inter-pocket 
impurity scattering but is still fairly visible for the values 
of impurity scattering used to fit NMR data. 
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